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ABSTRACT 

I discuss several topics in the applications of chiral symmetry at nonzero temper- 
ature. First, where does the rho go? The answer: up. The restoration of chiral 
symmetry at a temperature implies that the p and <li vector mesons are degener- 
ate in mass. In a gauged linear sigma model the p mass increases with temperature, 
VflpiT^ > mp(0). I conjecture that at the thermal p — tti peak is relatively 
high, at about ~ 1 GeV , with a width approximately that at zero temperature 
(up to standard kinematic factors). The OJ meson also increases in mass, nearly de- 
generate with the p, but its width grows dramatically with temperature, increasing 
to at least ~ 100 MeV by T^. I also stress how utterly remarkable the principle 
of vector meson dominance is, when viewed from the modern perspective of the 
renormalization group. Secondly, I discuss the possible appearance of disoriented 
chiral condensates from "quenched" heavy ion collisisons. It appears difficult to 
obtain large domains of disoriented chiral condensates in the standard two flavor 
model. This leads to the last topic, which is the phase diagram for QCD with 
three flavors, and its proximity to the chiral critical point. QCD may be very 
near this chiral critical point, and one might thereby generated large domains of 
disoriented chiral condensates. 

Based upon talks presented at the "Workshop on Finite Temperature QCD and 
Quark- Gluon Transport Theory", Wuhan, PRC, April, 1994 



1. Introduction 

In this paper I review several topics in which effective models are used to study 
the dynamics of chiral symmetry at nonzero temperature. The order is somewhat 
jumbled, and approximately in reverse chronological order from that in which the 
work was done. Whatever else, this has the virtue of presenting what I am most 
excited about (since it is most recent) first. 



2. Vector mesons and chiral symmetry 

The spectrum of dileptons in the collisions of heavy ions at ultrarelativistic ener- 
gies provides a window into "hot" regions of the coUison, whereby the formation of 
a quark-gluon plasma might be observed. For any effects of a quark-gluon plasma to 
be distinguishable from the background of ordinary hadronic processes, the system 



must last for a long period of time at a temperature T^, burning off the entropy of a 
quark-gluon phase. This is true for a first order transition with a large latent heat, 
but applies even if there is no true phase transition, as long as there is a large jump 
in the entropy in a narrow region of temperature. Numerical simulations appear to 
find such a jump in the entropy.El 

There are two distinct ways in which vector mesons at a temperature can be 
affected by the quark-gluon plasma, or more generally, by a hot hadronic plasma. We 
should speak of either type of plasma, since if there is no true phase transition the 
two cannot be distinguished, even in principle. (I speak only of nonzero temperature, 
and not of a system with nonzero baryon density. Cold, dense systems can be treated 
by an extension of the present methods, but may well involve new phenomenon.) 

The first is if the vector meson has a large decay width so that its lifetime is 
short, less than 1/m/c, and it decays inside the plasma. Then if its effective mass at 
Ty. is different from that at zero temperature, as it surely must be, in principle the 
shift in its mass, and so the corresponding peak in the dilepton spectra, is observable. 
At zero temperature the only example which appears in the dilepton spectra is the 
peak for the p meson, with a width of ~ 150 MeV. I suggested such a shift in the 
"thermal" p peak some time ago.! 

The second way in which a vector meson can be affected is if it has a long lifetime, 
much greater than 1/m/c, so that it decays outside of the plasma. Then even if its 
effective mass at is different from that at zero temperature, the overwhelming bulk 
of decays occurs outside of the plasma, and the peak for this vector meson in dileptons 
doesn't move. Instead, the peak shrinks, as the hot medium shakes the bound state 
apart. The J/\E' peak is suppressed in this way.i 

At first thought, one might expect that the u meson is like the J/\E' meson, 
since at zero temperature its width is small relative to that of the p meson, of order 
~ 10 MeV. In this work I argue that this expectation is wrong: the width of the u 
meson quickly grows with temperature, and is large at T^, at least ~ 100 MeV. Thus 
thermal u mesons decay inside, and not outside, the plasma, and the shift of their 
masses is in principle observable. The width of the p meson does not significantly 
increase with temperature. 

The first question to settle is: where do the p and u mesons go? Do their effective 
masses go up, or down, with increasing temperature? 

My first guess was "down" — that the effective mass decreased with temperature.! 
If the phase transition to a quark-gluon phase is primarily one of deconfinement, then 
this may be modeled by an effective bag constant which decreases with temperature. 
The effective mass of the p meson, like all hadronic bound states, should then decrease 
with temperature .H'ffl 

Numerical simulations of lattice gauge theory, however, demonstrate that the 
phase diagram of QCD with three flavors of quarks is rather involved, and depends 
crucially on the values of the quark masses. a This material will be discussed at greater 
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Figure 1: The phase diagram of a gauge theory with three colors and 2 + 1 flavors, 
following the results of the Columbia group. The quark masses are held in flxed ratio, 
so up = down quark mass, A4, and the strange quark mass = 32 A4. The |/-axis is A^, 
running from to oo, while the x-axis is the temperature T. There are two critical 
points, the deconflning critical point, Aio, and the chiral critical point, Aic- As 
drawn, QCD appears to lie above but near the chiral critical point. 

length in sec. 3, so I just summarize the results here. For simplicity I work with 
"2+1" flavors, holding the up, down, and strange quarks in flxed ratio, M. = Aiup = 
■M-down = rM. strange-, with T ~ 1/32.1 At iuflnitc M. there is a deconflning phase 
transition of flrst order, which might be — but at least to my eyes, today, need not 
be — modeled by a decreasing bag constant. 

(Perhaps it is worth elaborating on my concerns. At inflnite M. we are dealing 
with the phase transition in the pure glue theory. This is flrst order, with masses 
which shift with temperature on either side of the phase transition. But why should 
the masses decrease as the transition is approached from below? In particular, I 
can certainly rule out any transition in which the effective bag constant vanishes 
identically. If this were to happen, then by construction an infinite number of bound 
states would become massless. But there is no known critical point where an inflnite 
number of states become massless: for all known critical points in 3 + 1 dimensions, 
only a flnite number of flelds become massless.) 

Going down from inflnite M. there is a line of flrst order deconflning phase tran- 
sitions; for three colors and 2 + 1 flavors, this line of flrst order deconflning phase 




transitions does not extend indefinitely, but terminates in a deconfining critical point, 
Mv 

The opposite limit of = is the point of chiral symmetry for three flavors; here 
the chiral phase transition is expected to be of first order.Q Again, as Ai increases 
from zero there is a line of first order chiral transitions, which for three colors and 2 + 1 
flavors terminates in a chiral critical point, Aic- Thus the crucial feature of the phase 
diagram is a gap, for Aiv > Ml > Mlc, in which there is no true phase transition. 
Current lattice dataBl indicates that while QCD lies in this gap, Mv > Mqcd > Mc, 
Mqcd lies much closer to the chiral critical point than to the deconfining critical 
point. Thus QCD appears to be described by a smooth transition which is dominated 
by the (approximate) restoration of chiral symmetry. This phase diagram is illustrated 
in fig. (1), and will be considered at length later. 

2.1. Chiral symmetries 

Consequently we require a model which treats the restoration of chiral symme- 
try for both spin zero and spin one mesons. Before constructing the model, it is 
worth reviewing how the underlying quark fields, and so the composite meson fields 
constructed from them, transform under the global chiral symmetries. 

For simplicity I only consider the case of two flavors; since the meson is known 
experimentally to be predominantly ss, this is probably adequate for treating p and 
uj mesons. It is surely inadequate for mesons such as the r] and r]', which are not 
flavor, but SU{3), eigenstates. 

For left and right handed quark fields qi and qr, qi^r = (1 ±75)q'/2, under a global 
chiral symmetry transformation in the group SUi{2) x SUr{2) x f/a(l), 

qi e''" Ui qi , e'^^" Ur . (1) 

Here U^r are transformations under SUi^r{2), and e*^" that for axial f/a(l). 

The spin zero mesons constructed from the quark fields include two singlet fields, 

= qq , rj = q-f^q , (2) 

and two isotriplet fields, 

cr^ = qtq , vf = q-f^tq , (3) 

where the t are proportional to the Pauli matrices for SU{2). The rj and vr's are 
= 0~ states, while the and the cr^ are = 0^. The identification of the 
0"*" states is highly controversial; the candidates in the Particle Data Tables are the 
/o(975) and the ao(980). Both of these states have narrow widths, each less than 
60 MeV, and so are very possibly some other type of state, such as a KK molecule. 



The spin one mesons are constructed analogously by adding the Dirac matrix 7^ 
everywhere. In terms of the plausible candidates in the Particle Data tables these 
fields are the isosinglet states, 

= , /r = qi^rq. (4) 

and the isotriplet states, 

pi' = qYtq , af = qi^l'^tq- (5) 

The ti;(783), the p(770), and the ai(1260) are all familiar; the only unfamiliar state is 
the /i, which I identify with the lightest fi in the Particle Data Table, the /i(1285). 
The uj and p are = 1~; the /i and the ai are = 1+. 

I would like to make a side remark. From a study of the QCD phase diagram as 
a function of the current quark mass A^, it is known that how far QCD is from the 
chiral critical point depends crucially on the value of the mass for the isosinglet O"*", 
the (T^, at zero temperature.l Of course in QCD this is inevitably complicated by the 
decay channel into two pions, by mixing with KK molecules, and a myriad of other 
details. 

But in the quenched approximation, all of these details are absent, since by con- 
struction all mesons have zero width. Thus one can ask: 

In the quenched approximation, is the mass of the cr^ less than, or greater than, 
the mass of the p ? 

I am not suggesting that the answer in the quenched approximation is the same 
as for QCD; it could well differ. I am suggesting, however, that it is a well defined 
question, with an answer amenable by present day techniques. 

To return to the question of chiral symmetry, under SUi{2) x SUr{2) transforma- 
tions the scalars mix with each other as: 

an ^ n , 7] ^ , (6) 

while the vectors mix as: 

& /i, invariant ; p ^ ai . (7) 

For the scalars, the 0"*" and the 0~ states mix, while for the vectors, while the 1~ 
p and the 1"*" Oi mix, the singlet states, the u and the /i, are left unchanged. The 
singlets are invariant because they correspond to currents for fermion number and 
axial fermion number, which are conserved classically in the massless theory. From 
Eqs. (|^) and (J^), in a chirally symmetric phase the masses of the and the vr, the rj 
and the a^, and the p and the ai, are equal to each other; there is no prediction for 
the masses of the u and the /i, since they don't mix with any other states under the 
chiral symmetry. 



The axial Ua{l) symmetry is broken quantum mechanically. Nevertheless, for 
completeness I list how the states transform under Ua{l)'- 

ar, ^ 1] , ^ TT . (8) 

In contrast to the scalars, all vectors — the uj, /i, p, and the ai, are invariant 
under Ua{l) transformations. This follows because of the conservation of Ua{l) at the 
classical level. 

At zero temperature we know that the effects of the anomaly, and so the quantum 
mechanical breaking of the axial Ua{l) symmetry, are large, because the rj' (which 
for two flavors is the r] meson) is heavy. It is reasonable to suspect that at very 
high temperature the axial Ua{l) symmetry is at least partially restored. One way of 
seeing this is to see if the masses of the states which mix with each other under Eq. 
d^) are approximately equal. While this is possible for the singlet states, the r] and 
the arf, it is probably much easier for the isotriplet states, the vf and the cTt?. 

Thus on the lattice, one can indirectly look for the restoration of axial Ua{l) by 
measuring the mass splitting between the vr 's and the a-^ 's. This provides another 
motivation for measuring the properties of the scalar particles. 

In this paper I will ignore the possible effects of the restoration of axial Ua{l), 
and concentrate on the restoration of SUi{2) x SUr{2). For the scalars, this implies 
that I forget the rj and cr,r fields, keeping only the n and o"^ mesons. Since the there 
is just one sigma field, the isosinglet a^, I will refer to that simply as the a field. I 
stress, however, that even for two flavors there is another isotriplet of scalar fields 
which in principle should be included, the a^f's. As will be discussed in sec. 3, for 
three flavors there is a full nonet of sigma mesons which go along with the usual nonet 
of pseudoscalar Goldstone bosons. 

If we only wanted to decsribe the lightest fields at zero temperature, which are the 
pions, then it would suffice to use a nonlinear sigma model. While this is currently 
a very popular approach, I stress that at best it is most awkward to describe the 
phase transition to a chirally symmetric phase. The reason is rather trivial: above 
the temperature for the restoration of chiral symmetry, all particles must fall into 
degenerate multiplets of SUi{2) x SUr{2), and possibly Ua{l) as well. It is simply 
easier if we start with a theory in which all the fields which we need to form the 
complete multiplets are there to begin with, rather than having to generate them 
dynamically. Of course if we don't put the full multiplet in by hand, they will be 
generated dynamically. For example, in the nonlinear sigma model in 2+e dimensions, 
by using the large N expansion it can be shown that in the chirally symmetry phase 
the a field is generated dynamically, as a bound state of two pions.i 

Part of the prejudice against the hnear sigma model, as opposed to the nonlinear 
model, is due to the fact that the sigma meson is very broad at zero temperature, from 
its decay mode into two pions. Thus the feeling is, why bother? But even if the sigma 
is broad at zero temperature, inevitably it will become narrow at high temperature. 
This is trivial in a chirally symmetric phase, since then the sigma and the pions are 



degenerate. When the pion is massive at zero temperature, a narrow sigma can even 
emerge before the temperature of (approximate) chiral symmetry restoration. This 
happens because as the temperature is raised, the sigma mass goes down, while the 
pion mass goes up; thus the (thermal) sigma mass falls below twice the (thermal) 
pion mass before the temperature of chiral symmetry restoration.^ 

Thus the sigma meson(s), which is a broad state, lost in the hadronic mud at zero 
temperature, eventually emerges as a clean, nonzero temperature. 



2.2. Sigma models and vector dominance 



Having fixed upon a linear sigma model, we then have to include the vector 
mesons. Here we can rely upon one of th£_^umustly forgotten triumphs of the 1960 's, 
which is the model of vector dominance ;li2ryj indeed, it was Sakurai's gauge theory 
of the p which lead to the standard model, and thence to QCD. I hope to empha- 
size, however, that whatever its historical antecedents, that vector dominance is an 
amazing thing indeed. 

From the SU{2) Pauli matrices a"", I introduce t°- = cr"/2 and t° = 1/2; these 
matrices are normalized so that tr{t"-t^) = 6°'''/2, etc. Since we are forgetting about 
the rj and the a^^ fields, for the scalar fields we only need an SU{2) x SU{2) ~ 0(4) 
vector, 

$ = (Tt° + ivf-f. (9) 
Under a SUi{2) x SUr{2) chiral rotation, $ transforms as 

$ ^ nj^ilr . (10) 



For the vector fields, I introduce left and right handed fields as 

= (uj^ + /f) + (p^ + af) ■ r, 

= {uj^-m + (p^ - an-r. (n) 

The central question is how the vector fields transform under chiral rotations. The 
obvious, and certainly the most natural guess, is that they transform as one would 
expect under a global chiral rotation, which is homogeneously: 

Ai ^ nl Ai ni,r . (12) 



I now show that while the transformation of Eq. ( [T^ ) is standard, it does not lead to a 
model with vector dominance; the following repeats the discussion of Gasiorowicz and 
Geffen.ty I introduce the abelian field strength tensor for the left and right handed 
fields, 

J^l r \not VDM = 

d'Al - d'^Al , (13) 



and introduce the effective lagrangian for the gauge fields, 



^notVDM — -tr y{Fl^'^\notVDM) + (F'r'^lnotVDM) ) 

+ tr [{Aff + iA>^f) (14) 

This is a reasonable effective lagrangian for the vector fields. It respects the global 
chiral symmetry of SUi{2) x SUr{2), while the mass term for the gauge fields will 
ensure that the gauge fields are massive. Of course we should also add the coupling 
to the scalar field $, but let us neglect that for the time being. 

The current for Cnot vdm is computed by the standard Noether construction. For 
the left handed fields, say, take f2/ = exp{iui-i^, and compute the infintesimal variation 
of the lagrangian for a spatially dependent uf, then the current is the coefficient 
of d^uji. For infintesimal ui, in group space the change in the vector potential is 
6Af = [Af, ui]. Immediately we see that the mass term in Cnot vdm doesn't contribute 
to the current, since 6tr{{Afy) = tr{Af[Af,u!i]) = by the cychc property of the 
trace. That doesn't mean that there isn't a current; the abelian field strength tensor 
gives rise to a perfectly fine (left handed) current, 

flnotVDM = [Ald^A\-d^A^]. (15) 

From the standpoint of a current, this is a very reasonable expression, bilinear in the 
fields. This is analogous to the current we know for a scalar field, 0*9^0 — (9^0*0, and 
that for a fermion field, ip'f^ip. 

But it's not the current required for vector dominance. Vector dominance is the 
statement that the largest coupling of hadrons to photons is through a current linear 
in the appropriate vector fields, with a coupling proportional to the mass (squared) 
of the vector fields. We are not used to dealing with currents linear in the fields! 

Such a current can be constructed, but at the expense of quite a leap. After 
all, we are dealing with a theory which has only a global chiral symmetry. Vector 
dominance instructs us to promote this symmetry to one which is local. At least for 
myself this was a remarkable step which I strenously resisted, to no avail. The basic 
point can be understand by changing the transformation of Eq. (|T^) to that for a 
local chiral rotation, 

K - K ^l,r . (16) 

-ig 

1 introduce the covariant derivatives Dj^^ = — W^ir ^^"^ coupling constant g. 
Notice that the coupling constant g is that for vector dominance — it has nothing 
(directly) to do with the coupling constant of QCD, and is generically a large coupling, 
of order one. 

The amazing point is that by assuming that the vector fields transform mhomogeneously 
under (local) chiral rotations, we are automatically guaranteed that the mass term 



gives the proper current. This is simply because for small ui, SA^I = i d^uji/g+ [AJ*, uji\ . 
Now for the mass term, the commutator term doesn't contribute to the current, as 
before. But the inhomogeneous part of the transformation of the gauge field then 
gives, trivially, a contribution to the current as 

2 

Tfl 

Jt = y • (17) 

This is the type of expression required for vector dominance: the current is linear in 
the gauge fields, with a coefficient proportional to the mass (squared) divided by the 
vector coupling constant. 

Since this is the current we want, we then have to ensure that no other terms 
in the lagrangian contribute to the current. But having taken the giant step of 
introducing a local chiral symmetry, the rest is easy. To ensure that no other terms 
in the effective lagrangian contribute to the current, we require that all other terms 
couple to the vector fields in a gauge invariant manner. As long as the couplings 
are gauge invariant, then, by definition the rest of the lagrangian will be invariant 
under the gauge transformation of Eq. (p!6D! Hence we introduce the nonabelian field 
strengths, 

= d^Al - d'^Al - ig [Al, A^] . (18) 
For the scalar field $ the appropriate covariant derivative acts as 

= d^"^ - ig (Af<l> - $A^) . (19) 

There is a relative minus sign between the coupling constant for the left and right 
handed handed fields because from Eq. ([T0|) , $ transforms under fir and Q]. Given 
these quantities, the effective lagrangian for the gauged linear sigma model which 
respects vector dominance is then 

C = tr\Df,^\'^ - htr{^) + tr\^\^ + ^ {tr\^\'^f 

+ \tr {{Frf + {Frf) + rn' tr [{A^f + (Aff)^) (20) 

The scalar potential is standard, with a term linear in $, ~ htr[^), added to ensure 
that the pions are massive. 

The assumption that the spin one mesons couple to the scalar field only through 
coupling which are locally gauge invariant greatly restricts the possible couplings. If 
only a global chiral symmetry is imposed, many more terms are possible. For example, 
the terms $ and ^A!^ both transform under global chiral rotations like $ itself, Eq. 
(p!0|). Since QCD preserves parity, we can require that all possible terms be invariant 
under the interchange of left and right handed gauge fields. For example, for the 
quartic terms, besides |y4f$ - $y4^p, the terms + and |v4f$p + 

are also allowed under the global chiral symmetry. 



From the viewpoint of the renormahzation group, this restriction in possible terms 
is most striking. The standard prescription of the renormahzation group is the fol- 
lowing. Start with the terms with the lowest mass dimension, which are usually the 
mass terms. Write down all mass terms consonant the relevant symmetries. For 
terms with higher mass dimension, which are typically cubic or quartic in the fields, 
and whose coupling constants are less relevant, all possible couplings which respect 
the symmetries are allowed. For example, for a theory with an N component vector 
0, if the symmetry is 0{N), only terms as ■ arise. But if the mass term is only 
invariant under 0{N — 2), many more terms are possible, as long as they are 0{N — 2) 
invariant. 

What is so peculiar about the model of vector dominance is that the mass term 
for the gauge fields is the most relevant operator and manifestly breaks the gauge 
symmetry. Yet for all other terms — the coupling of the gauge fields to themselves, 
and the coupling of the gauge fields to the scalar fields — one uses this local gauge 
symmetry, most crucially, to greatly restrict the possible couplings of the gauge field. 
Of course what is unusual about the model of vector dominance is that it is a local, 
and not simply a global, symmetry which arises. Somehow this must be crucial to 
its use and eventual consistency. But it seems to me as if something more interesting 
than mere dusty phenomenology is afoot. 

1 only use the effective lagrangian of Eq. (^) at the tree level. One simplifica- 
tion follows immediately: the equations of motion for the gauge field are DiFi"^ = 
m'^Af + ji, where jf is the current from the scalar field. To lowest order this re- 
duces to d^Fj^^ = m'^A'^, which is only consistent is the gauge field is in Landau 
gauge, d^Ai = 0. Thus in perturbation theory the propagator for the gauge field is 
[6'^'^ — p'^p'^ /p'^)/{p'^ + m^), which falls off like ~ at large momenta. Of course in 
gauges other than Landau, the propagator behaves like {6'^'^ — p^p^ /m?)/ {p^ + m^), 
and is of order ~ 1 at large momenta. So at least the ultraviolet behavior of the 
gauge propagator isn't as horrible as it might be. 

Of course there is still the question of how to consistently quantize the theory. 
The usual linear sigma model at least has the virtue of perturbative renormalizability, 
although one may well be pushing it into a highly nonperturbative regime. The 
effective lagrangian of Eq. (^OD does not have such a virtue, however, and it is far 
from clear how to treat the standard problems of such theories: the coupling of 
(unphysical) massless modes which are present in the gauge propagator to physical 
modes, unitarity, lack of renormalizability at two loop order, and so on. It may well 
be possible to find some extension of Eq. (pO]) which is (at least) renormalizable, but 
I haven't thought seriously about it. 

Much of the physics of vector dominance can be read off from the coupling of the 
spin one fields to the scalars through the covariant derivative. After working out the 
matrix algebra, one finds 



+ I (^'^vf + g{p^ X Tf) - gaarf + C + fi ■ (21) 



Notice that the a; meson drops out completely: from Eq. (|llD, in the vector fields Ai^r 
the u meson is proportial to the unity matrix, and so the u cancels in the covariant 
derivative in Eq. (pj]). There is a coupling of the cu meson, due to the anomaly, but 
I will only discuss that in passing here. 

The principal couplings of interest in Eq. ( pTj ) are the coupling of the p to two 
pions, as p'^ ■ (vr x S^vf), and the coupling of the ai to the a and vf, a a^^ ■ vf. Both 
of these couplings are proportional to the coupling of vector dominance, g, and are 
presumably responsible for the principal decays of these particles. 

To obtain more realistic expressions, I assume that the potential of the scalars 
is such that a vacuum expectation value for the sigma field develops, p^ = p^ < 0. 
Because of the background magnetic field ~ h, the vacuum expectation value of 
must be along the a direction. After shifting a ^ + a , however, a complication 
arises which is special to the gauged sigma model: there is a cross term between the 
ai and the pion, as ^af ■ d^ii. This cross term must be eliminated by a shift in the 
ai field, 

9(^0 gf,^ ^ ^22) 



One can show that the equations of motion for the shifted ai field are such that the 
shifted (and not the unshifted) ai field should be in Landau gauge. 



After this shift, for the vector fields the mass term written in Eq. (|20|) , tr{A'j 
gives 

ml = ml = , m^^ = mj^ = + (gao)^ . (23) 

This is in good accord with experiment, where the p and the u, as well as the ai 
and the f\, are very nearly degenerate. Yet I find it mystifying. Besides the mass 
term tr{Afj.), there is a second set of mass terms which are invariant under global 
SUi{2) X SUr{2) rotations, {tr{Ai^r)Y ■ These new mass terms only contribute to the 
masses of the u) and the /i, but experimentally they do not seem to be there: the 
masses of the u and the p, and the masses of the ai and the /i, are each within 2% 
of each other! Perhaps my understanding of vector dominance is faulty, but there 
appears to be a further principle at work, something that tells us that such flavor 
singlet mass terms are very small. Of course such singlet mesons can only be treated 
realistically within the context of a full, three flavor model, but I believe that the 
paradox will remain. 

After shifting the ai field as in Eq. (0), the kinetic term for the vr becomes 

^ ^ ^ (5M^)2 _ (24) 



2 {m"^ + {gaof 



This demonstrates what is known as the "partial" Higgs effect. In the limit that 
the mass of the vector meson becomes very heavy, m — > oo, all effects of the vector 



mesons should decouple, and the pions (sigmas, etc.) are unaffected. The opposite 
limit, when the mass of the vector meson vanishes, m — > 0, is also familiar. From Eq. 
(p^), the kinetic term for the pions vanishes as m — > 0, but this is simply because 
in this limit there is a true gauge symmetry and a true Higgs effect, with the pions 
turning into the longitudinal components of the Oi's. For finite but nonzero m, what 
happens is that ratio y^m^ + [gaoY/m = ma^/nip enters. Experimentally this is a 
significant number, ~ 1.6. 

In this vein, I should also mention what is known as the KSRF relation.0 This 
relation predicts that the ratio of the Oi to p mass is fixed, = \/2. Within the general 
context of gauged linear sigma models, there is no reason why the KSRF relation 
should be satisfied; see, for example, the comments of Lee and Nieh.li2l One could 
argue that the actual value of ma^/mp ~ 1.6 is not too far from \/2 ~ 1.414..., but I 
see no particular virtue in adhering to the KSRF relation. Instead, I use the ma^/mp 
ratio to fix the value of the vector meson coupling g. 

After rescaling tt — > {maj^/mp)'n', the physical quantities of interest are: 

ml = - , ml = ^ + 2\al , /. = ^ a, . (25) 

\mp ) (To (To TTLa^ 

Fitting the above expressions with = 137 Me V", mo- = 600 Mel^, rup = 770 Me V", 
and rriai = 1260 MeV, the results are 

ao = 152 MeV , h = (102 MeV f , 

= -{588MeVY , A = 7.6 . (26) 

For comparison, if we send the mass of the vector mesons to infinity, and denote the 
analogous quantities by a superscript nv for "no vectors" , with the same masses for 
the TT and the a I find 

a^"" = 93 MeV , h""" = {120 MeV f , 

(/i™)2 = -{567 MeV f , A™ = 19.7. (27) 

From Eq. (|25|) , the "partial" Higgs effect tends to increase the pion mass and 
to decrease the pion decay constant f-,,. This description is somewhat misleading, 
because it is phrased in terms of the quantities h, ao, and /Xq, which are not directly 
physical; the physical quantities are the masses and /t^. One quantity which is physical 
is the value of the scalar self coupling. A; that is a dimensionless number, which tells 
one how close the linear model is to the nonlinear model, the nonlinear model being 
recovered as A — oo. I find it notable that adding the vectors changes this coupling 
from A"^ = 19.7 to a much smaller value, A = 7.6. Presumably with the p and ai 
vector mesons, some of the scalar interactions, which otherwise arise solely from their 
self interaction, A(|<l>p)^, can be made up by exchange of vector mesons, giving a 
smaller value of A. 



In conclude this section by discussing how vector meson dominance works in the 
hmit of a large number of colors, A^^^ oo; if the nonabelian coupling of the SU{Nc) 
gauge theory is (7^, g'^Nc = is held fixed and of order one as ooM Consider, in 

particular, the spectrum of dileptons, which is given by the imaginary part of the two 
point function of hadronic currents, {j'^{x)j^{0)) . In perturbation theory this graph 
starts with a free quark loop, which is of order Nc for quarks in the fundamental 
representation of SU{Nc), plus an infinite series of corrections in as- Now at large 
Nc it is known that all color singlet currents are saturated by hadronic states, where 
mesonic masses are of order one. This works fine for the vector channel considered 
here if the coupling of to the p is of order \/Nc, which is standard at large Nc- 

This can also be seen from vector dominance, assuming that vector dominance 
holds in the large limit. From Eq. (|17D, = im? A^/g\ m = rrip is of order one, 
as is the vector field A^^, since it is just another mesonic field. The coupling of vector 
meson dominance, however, is a trilinear coupling between three mesons. This is 
typically vanishes at large Nc as g l/\fWc, so Eq. ([T7| ) becomes ~ ^/WcA^, and 
agrees with the counting of quark loops. 

All of this is totally standard. Let us then push on and consider corrections to 
the large Nc limit. The coupling of the p to two pions is again a trilinear meson 
coupling, of order g ~ 1/a/]Vc, so the width of the p meson is of order g"^ ~ '^/Nc- 
Thus while the height of the p peak is very large, of order Nc-, it is very narrow, of 
order l/Nc-, so the total area under the p peak is of order one. At order one there 
are other contributions to dilepton production, such as from pion pairs. Since the 
coupling of hadrons to photons is electromagnetic, and pions don't interact at large 
Nci the pion contribution to dileptons is just given by the imaginary part of the one 
loop diagram of pions. Since both the electromagnetic coupling of the pions and their 
number are of order one, so is the total pion contribution. Thus the large Nc limit is 
consistent with vector meson dominance as follows: vector mesons, such as the p, 
etc. give tall but narrow peaks, under which is a smooth continuum for 7r"^7r~, K^K~ 
pairs, and so on. The total area under the vector meson peaks, and the continuum 
from scalar meson pairs, is the same, of order one. 

The above counting just shows that the large Nc limit is compatible with vector 
meson dominance; it does not show that in fact vector meson dominance emerges 
at large Nc. Nevertheless, it is important to see that there is nothing inconsistent 
between vector meson domainance and the limit of large Nc. 

2.3. Vector mesons at nonzero temperature 

With all of this introduction, the extension to nonzero temperature is relatively 
straightforward. I deal exclusively with mean field theory, since at least that gives 
a well defined approximation. I am not claiming that it is a good approximation, 
simply that it is well defined. 



Consider the mass terms in the effective lagrangian of Eq. (pO|), 

C^ass = /i'tr|<l>|2 + mHr({A'if + {A'^f) (28) 

Diagramatically, both masses increase with increasing temperature. For example, from 
the scalar self interactions, a term as ~ AT^tr|$p is generated from the interactions 
with a pion loop; similarly, the interactions of vector mesons with a pion loop also 
produce a term ~ g'^T'^tr{Af + Af.). A priori, I do not see any reason why the 
contribution of a pion loop to one term should be significantly larger, or smaller, 
than the other. 

This allows me to make what is mathematically trivial, but physically profound, 
observation: that the effects of thermal fluctuations are always to cause both /i^ and 

to increase with increasing temperature. The details of their relative increase 
surely depends on details of the model — the number of flavors, colors, the values of 
the self couplings, etc. But in any gauged sigma model, both mass terms go one way 
and only one way with increasing temperature, and that is up. 

The expressions for the masses of the vector mesons given in Eq. (|23|) remain 
valid as and /i^ (and so cTq) change. Thus a gauged sigma model predicts that the 
mass of the p increases monotonically with temperature; similarly, that the masses 
of the u! and the also increase. Because of the constraints of chiral symmetry, the 
masses of the other vector mesons, the ai and the /i (and their strange partners) 
behave in a more complicated fashion: while increases, /i^ becomes less negative, 
so (To decreases. Thus the m^^ and m/^ first decrease with temperature, until they 
become (approximately) degenerate with the p — uj; then the masses of all states 
continue to increase with temperature. 

The above assumes, implicitly, that there are no singlet masses for the vector 
fields, (tr(y4i j,))^. Such singlet mass will split the uj and /i apart from the p and ai. 
Even if such a term is not present at zero temperature — because of experimental 
constraints on the known masses of the p and u — unless there is a symmetry reason 
why such a term cannot arise, it will be generated at nonzero temperature. If such 
a term is generated, the degeneracy between the u and the p will be lifted, but both 
masses will still increase with temperature. In any case, as we shall see, all states 
becomes broad at nonzero temperature, so a relatively small shift in the real part of 
the masses may not be significant. 

The prediction that the mass of the p (and those of the u and (p) increase with 
temperature is not standard lore. It is a prediction of a gauged sigma model, al- 
though the generality of the remark does not seem to have been stressed in previous 
studies.0 Different models, such as those of Brown and predict the opposite: 

that the mass of the p decreases as the temperature goes up. In principle this could 
be studied by means of other phenomenological approaches, such as sum rules. Here 
one runs into a problem: extending sum rules to nonzero temperature is not free 
of ambiguity. Thus while some papers find that the mass of the p decreases with 



increasing temperature others find that it increases. Besides my obvious preju- 
dice in the matter, it is worth noting that the studies which claim to find that the 
p mass increase with temperature find that it does so in a manner consistent with 
the restoration of chiral symmetry. Of course chiral symmetry alone is not enough to 
claim that the p mass goes up with temperature; all chiral symmetry tells you is that 
the p and oi should be approximately degenerate. 

One approach which differs from a gauged linear sigma model, and yet incorpo- 
rates chiral symmetry in a crucial way, is the study of Nambu-Jona-Lasino (NJL) 
models.0 If the predictions of the gauged linear sigma model have any generality 
whatsoever, they should hold for a NJL. There is a technical caveat: what is usually 
studied are (chirally) symmetric NJL models with the simplest possible interaction, 
{tptil^Y + (ipt'y^ip)^. This type of interaction is certainly adequate for treating the 
scalar and pseudoscalar particles, since undoing the quartic interactions by means 
of auxiliary fields with automatically generate effective fields for these particles. To 
study vector and axial vector particles, however, I suggest that it is crucial to include 
interactions such as [ipt'^^ipY + i what is termed the extended NJL 

model. 

With this cautionary aside, there is a unique prediction: in any extended NJL 
model, the mass of the p should increase, monotonically, with temperature. If not, 
there is something seriously wrong with a gauged linear sigma model, or with vector 
meson dominance. 

In principle, the masses of the vector mesons could also be studied by numerical 
simulations of lattice gauge theory. This is not elementary to date all studies have 
concentrated on the (manageable) task of computing static correlation lengths. While 
the TT and a channels are special, in all other channels one finds that two quark states 
have a static correlation length ~ l/(27rT), three quark states one of ~ l/(37rT), 
and so on. The usual interpretation is that this represents the propagation of free 
quarks, dressed by interactions. I suggest that if the deconfining transition is far 
afield, then what one is seeing really is the propagation of mesons (and baryons) in a 
hot hadronic gas. To definitively answer the problem, however, it is necessary to look 
not at static correlation lengths, but at the true poles in the appropriate (effective) 
propagators. These are defined in real time, by the analytic continuation of the 
spectral densities from euclidean momenta. Thus it is necessary to measure not just 
the static correlation lengths, at p° = 0, but the correlation lengths for p° = ±27rT, 
±47rT, etc., and then fourier transform. This is a very difficult problem. 

At this point I note an alternate scenario by Georgi.El He works in a nonlinear 
model, with the term for the p mass squared proportional to a dimensionless coupling 
constant times He then argues that as 0, nip 0. In terms of the hnear 

sigma model, it is not apparent to me why the mass term for the p should involve /^r 
at all; why is it simply not some other dimensional parameter? One can rephrase the 
argument in terms of the linear sigma model: Georgi's term corresponds to assuming 
that there is no bare mass term for the p, such as m?tr{A^ + A^), but only a term 



as tr{^^^)tr{Af + Af.), which would have a dimensionless couphng constant. In the 
hnear sigma model this is odd: why should there be no bare mass term for the vector 
mesons, yet one induced by spontaneous symmetry breaking? After all, the bare 
mass term is more relevant than Georgi's term. Leaving prejudice aside, I appeal 
to the lattice data. Admittedly, while static screening lengths are not automatically 
the (inverse) pole masses, it is true that if the p were massless at T^, then its static 
screening length would diverge, like that for the a and the vr. This is not seen in 
present day lattice simulations: only the a and the vr have divergent static screening 
lengths (at least for two flavors, where the transition appears to be of second order). 

Returning to the matter at hand, one of the drawbacks of the gauged sigma model 
is that it does not predict how much fi^ and increase with temperature. I thus 
adopt what is, most honestly, a wild guess. As the temperature increases, the mass of 
the p goes up, while that of the ai goes down, until they meet (approximately), after 
which they both increase. I define the temperature of chiral symmetry restoration 

as that for which the two masses are approximately equal, and assume that they 
obey the rule that their mass at is the arithmetic mean of their masses at zero 
temperature, 

mp(T^) ^ m„,(T^) ~ ^ (mp(0) + m,,(0)) = IGeV . (29) 

I emphasize that this is only a guess; at least it is easy to remember. It probably 
is an upper bound on how much the p mass can increase with temperature. For 
example, Son^l^ has done calculations in a gauged nonlinear sigma model, and finds 
that with different parametrizations, the mass of the p can even stay (almost) constant 



with temperature, or increase as much as in Eq. (^9]). The former seems extremely 
unlikely to me; again, the pion loops that drive /i^ up also contribute to driving 
up. Nevertheless, how much the p mass increases with temperature is probably only 
something which can be settled by lattice simulations. 

The virtue of Eq. ( P^ is that it gives a unique parametrization of how masses 
change with temperature; I find that if pi = — (588Me\^)^ and ml = {770 MeV)'^ 



are the values at zero temperature, then Eq. ( P9D is satisfied for 

= rnl + .6 ill^^ . (30) 

One can then compute the properties of the model as a function of a single parameter, 
p"^] all effects of increasing temperature are then modeled by increasing /i^. Again, I 
stress that Eqs. ( PU[ ) and ( |5(]| ) are only guesses. 

In fig. (2) I give the values of the masses of the vr, a, p, and ai as a function 
of increasing /i^ (~ temperature squared). There is no prediction whatsoever for T, 



X' 



since one can only see that there is first approximate chiral symmetry restoration for 
/i^ ~ +(150 MeK)^. AT T^, rria- ~ ~ 200 Me V^; by assumption, nip ~ ~ 
1 GeV. In fig. (3) I show versus /x^; the relevant observation is that at p'^, f-,, 
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Figure 2: Masses of the vector and scalar particles versus the mass parameter /x^ in a 
two flavor, gauged linear sigma model. The masses are plotted against the y-axis in 
MeV; from top to bottom, these are the masses of the ai, p, a, and vr, respectively, 
/i^ is along the x-axis: the values at zero temperature are at the extreme left, with 
increasing /i^ equivalent to increasing temperature squared. 

has decreased to about a third of its value at zero temperature. This decrease in f^^ 
is elementary: for vanishing pion mass /^r ~ o", and vanishes at the temperature for 
chiral restoration (if the transition is second order). 

We can use the effective model to make more detailed predictions about the 
widths of the states at T^. At present I give only a very crude estimate, with careful 
estimates given later.lll There is certainly new kinematics which opens up: at zero 
temperature I assume that m„ = 600 MeV, so the p can't decay into a's plus tt's. 
That changes at T^, since then rricr ~ "n^Tr, and furthermore the p is heavier. 

The decay modes of the thermal p can be read off from the effective lagrangian 
of Eq. (|2T|), but the allowed channels follow from general principles. Even if it is 
possible kinematically, the p doesn't decay into two cr's because of isospin symmetry; 
also, p doesn't decay into an because of G-parity. Of course at nonzero temperature 
p vrvr, as at zero temperature. By similar arguments one can see that the only 
new three body decay mode is p — rnia; this respects both isospin and G-parity, 
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Figure 3: The value of the pion decay constant, versus the mass parameter /x^ in 
a two flavor, gauged linear sigma model. The y-axis is the value of in MeV, while 
fi^ is along the x-axis: the values at zero temperature are at the extreme left, with 
increasing /x^ equivalent to increasing temperature squared. 



and with the assumed masses has nonzero phase space. However, I do not believe 
that this three body decay gives a significant contribution; the coupling constant 
for p — » vrvr is ~ ^; from Eq. (|2T|), that for p rnia is g times g"^ fT^rrifj /rn?p, or 
~ ^/3. Thus the amplitude for p —>■ vrvrcr is about 10% that that for p —>■ tttt, even 
without the further restrictions which arise kinematically from a three body, versus 
a two body, decay. For this reasons I neglect the three body decays of the thermal p. 
For the remaining mode of p — vrvr, I make life even easier for myself by neglecting 
thermal effects entirely, computing the decay as if it were at zero temperature. This 
is valid because the pions are so energetic; if n{E) = l/{exp{E/T) — 1) is the Bose- 
Einstein distribution function for the pions, at ~ 500 MeV and T ~ 150 MeV, 
2n{E/T) ~ .07. Consequently, 



r,(Tj ~ ^ r,(0) ~ 175 MeV. (31) 



As the mass of the n goes up, the decay width of the p goes down, to about 100 MeV 
for rup ~ 1 GeV and ~ 300 MeV. Given the crudeness of the model, I conclude 
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Figure 4: The width of the p, under the extreme assumptions discussed in the text. 
Essentially, the formula for the width is assumed to be the same as at zero temper- 
ature, with the only difference from the mass shifts of the p and tt. The y-axis is 
Vp{Ty?j in MeV , the x-axis the value of m.„{T^) in MeV. 

that the width of the thermal p at IGeV is approximately equal to its width at zero 
temperature. In fig. (4) I present a graph of rp(T^) for rup = IGeV, and different 
values of m^r, under the above assumptions. Due to phase space, the decay width 
decreases from that in Eq. (|3T|). 

Implicitly I am assuming that the coupling constant for vector meson dominance, 
g, does not change with temperature. This is probably wrong, but consideration of 
the one loop /^-functions suggests that g should decrease with temperature near a 
critical point, making the width even smaller. 

Computing the width of the u meson is more difficult. A proper analysis requires 
an understanding of the Wess-Zumino-Witten term,lli which I defer for another day.tl 
For now I make a very simple assumption, that the decay of the u is dominated 
by what is known as the Gell-Mann Sharp Wagner (GSW) mechanism. This is the 
statement that the decay ^ 3 vr proceeds hj u ^ pn, with the virtual p decaying as 
p — >■ TTTT. The amplitude for the first process, — > pvr, is dominated by the anomaly, 
and proportional to 1/ fn] the amplitude for the second process is proportional to g 
times kinematic factors. At presentil I ignore all of these details to concentrate simply 
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Figure 5: Masses of the scalar particles versus the mass parameter /i^ in a three flavor 
linear sigma model which is not gauged. The masses are plotted against the y-axis in 
MeV; /i^ is along the x-axis: the values at zero temperature are at the extreme left, 
with increasing /i^ equivalent to increasing temperature squared. Notice that unlike 
the two flavor model, for three flavors I take rricr^, (0) = 975 MeV. The assumptions 
which produce this plot are described in section (3.2.1). The point of note here is to 
see that the kaon increases in mass, to approximately 650 MeV. 

on the change in with temperature; since f-wiTy^) ~ /7r(0)/3 and r^(0) ~ 10 MeV, 



This estimate is very crude, and clearly requires a more careful analysis.li5l The basic 
point, however, is bound to be correct: as chiral symmetry is (approximately) re- 
stored, decreases, leading to a much broader uj meson than at zero temperature. 
Kinematic factors will only help, as the uj moves up in mass, so the allowed phase 
space increases as well. 

Lastly, although the model discussed does not incorporate it, it is reasonable 
to ask about the meson. The chiral partner of the 0(1020) is one of the two /i 
states, either the /i(1420) or the /i(1510). There is some controversy concerning the 
/i(1420), so I take the /i(1420) as some sort of special quark state, and assume that 




(32) 



the chiral partner of the cj) is the /i(1510). Using the same kind of simple minded 
prediction for the IHRSS Clt Ty^ as in Eq. (^), I suggest that 

m<^(Tj ^ nif,{T^) ~ ^ (m^(0) + my,(i5io)(0)) = 1250 Me^ . (33) 

Even with this generous increase in the mass of the thermal 0, the thermal K mass is 
like that of the vr, and also increases, to at least 650 MeV] see, for example, fig. (5).t3 
With these numbers, the does not have phase space to decay as (p ^ K K (except 
for thermal processes, which should be small). Even if the masses were to change so 
that (f) ^ K K were allowed, it is very unlikely for the width of the (j) to change so 
much that it decays inside the plasma. 

In summary, for the thermal vector mesons at the temperature of chiral symmetry 
restoration, T^, the width of the p remains about that at zero temperature, and decays 
inside the plasma; the uj becomes much broader, so much so that most thermal 
cj's decay inside the plasma; the remains narrow, decaying primarily outside of 
the plasma. I do not expect a more careful analysisli3 to significantly alter these 
conclusions. 



3. Disoriented Chiral Condensates 



If an extraordinarily large number of pions are produced in a hadronic collision, it 
is natural to think that these unusual events might proceed by the coherent decay of 
a (semi-) classical pion field.i2li There might be a very distinctive signature of such 
a coherent decay: a given classical pion field points in some direction (or directions) 
in isospin space, so coherent production is necessarily one in which isospin symmetry 
would appear to be, at least locally, badly violated. That is, most of the pions from a 
given region would be primarily charged, or conversely, primarily neutral. (Of course 
in the total collision, isospin will always, in the end, average outj Such behavior may 
have been observed in Centauro events in cosmic ray collisions.c3 Experimentally the 
situation is unclear: some, but not all, cosmic ray experiments see Centauro behavior, 
while collider experiments do not see Centauro type events. Hopefully the "Mini Max" 
experiment at FNAL may help settle the issue. In the interim, it behooves theorist 
to consider seriously the possibility. 

I introduce the quantity 7^3, which is the ratio of neutral to the total number 
of pions. In the average, hadronic collisions conserve isospin. This will be true if 
one averages over different events (say in pp collisions), or over all regions in rapidity 
of A A collisions. Thus in the average there will be a binomial distribution in TZ^, 
strongly peaked about the isospin symmetric value of 1/3. 

But the behavior of the average might not be representative of all events. Suppose 
that in individual events single domains are produced, in which the pion field has a 
nonzero vacuum expectation value which points in a given, fixed direction in isospin 
space. If all directions in isospin space are equally likely, while the average value 



of 7?.3 in that domain is 1/3, the distribution is far from binomial, = l/(2\/^).0 
At this point I haven't spoken of how, in detail, such domains could be produced 
dynamically; a possible mechanism is the subject of the following. 

Note that while I have chosen the 3 direction in isospin to be that for neutral 
pions, there is nothing special in this choice. For example, if TZi is the fraction of pions 
in the isospin 1 direction, then the distribution in TZi is l/{2^/Rl). This is obvious, 
since by an isospin rotation the 3 direction can be relabeled as the 1 direction. The 
detailed form of the distribution in TZ^ or TZi is a consequence of the symmetry group 
of isospin being 0(3), and differs for other symmetry groups. The important point is 
that the distribution is far from binomial. 

Bjorken, Kowalski, and Taylor have proposed a "Baked Alaska" scenario wherein 
such a distribution is produced in pp collisions; the idea here is to trigger specifically 
on pp collisions in which the multiplicity is very large. 011111 They refer to a domain 
with a fixed direction in isospin space as a "Disoriented Chiral Condensate" (DCC); it 
is disoriented because in the true vacuum the pion field doesn't point in this direction 
(instead the sigma field acquires a vacuum expectation value). In some abuse of 
terminology, I refer to any system which generates an isospin distribution =1/ {2\/lZ^) 
as that of a DCC. It is worth remembering, however, that the way in which DCC's 
can arise in hadronic collisions can, and in general will be, very different from that in 
AA collisions. 

The behavior of central AA collisions at high energies is manifestly an example 
of a system with a large number of pions; for example, at RHIC and LHC, there 
will be on the order of thousands of pions per unit rapidity. The hope is then that 
while the total isospin in conserved for the total collision, that for given slices in 
rapidity, the distribution ilL.'T^s is not binomial, but that of a DCC. As emphasized 
by Blaizot and Krzywicki,E2l the difficulty is understanding why domains should be 
large in the transverse dimension: many transverse domains, oriented in different 
directions, washes out the DCC effect to give the usual binomial distribution. 

Rajagopal and WilczekS (R&W) proposed a dynamical mechanism for obtaining 
DCC distributions in heavy ion collisions. They use a linear sigma model to describe 
the chiral behavior over large distances, and make a dramatic assumption about the 
dynamical evolution in time. They assume that the dynamics is that of a "quench": 
the initial state is chirally symmetric, as appropriate at high temperature, but its 
evolution forward in time is by means of of the (classical) equations of motions at 
zero temperature. R&W assert that this dynamics, which is very far from equilibrium, 
produces amplification and coherent oscillation of pions at low momentum; that is, 
large DCC domains. 

In this talk I review two attempts by S. Gavin, A. Gocksch and Lto understand 
how disoriented chiral condensates might arise in heavy ion collisions.c3Q'El Our work 
was directly motivated by the work of Bjorken, Kowalski, and Taylor, and also by 
the work of Rajagopal and Wilczek. At the outset I must confess that while we all 
speak of DCC's, in detail the mechansims of Bjorken and of Rajagopal and Wilczek 



are really very different. In the first section we consider the work of Rajagopal and 
Wilczek in detail, extending their results.l^ Unfortunately, our results are negative: 
we find that the only way to get large DCC's is if there is at least one light field 
about. In the third section I turn to what seems a long digression on a disconnected 
topic: the phase diagram of QCD with 2 + 1 flavors. This leads, however, to the 
speculation that the equilibrium phase diagram for 2 + 1 flavors could itself generate 
a large distance scale. This then might generate large domains of DCC's. 



3.1. DCC's in a two flavor model 

In this subsection I review resultsii which extend previous work by Rajagopal and 
Wilczek, R&W.EZI Our major purpose was to concentrate on the question of domain 
size and on the experimentally relevant question of the nature of the distribution in 
7^3. While our results, like those of R&W, are the product of numerical simulations, 
much of our understanding follows from the analysis of Boyanovsky et al.f^ who 
consider the question of domain growth following a quench. While the work of (30) 
was done primarily in weak coupling (as applies in inflationary cosmology), it is direct 
to extend it to strong coupling, at least qualitatively. 

For the purposes of discussion we take the assumptions of R&W for granted. 
This means that we consider only the effects of two quark flavors by means of a linear 
sigma model, ignoring vector mesons. This involves an 0(4) vector $ = (cr, vf) where 
a is an isosinglet = 0^ field. The action for the $ field is 

£ = 1 (9,$)^ -//.$ + ^2^2 ^ A^^2)2_ 

At zero temperature the vacuum state is ($) = {v, 0). The parameters of the model 
are the quartic coupling A, the vacuum expectation value of the sigma field, v, and 
a background magnetic field, H = h{l,0), which makes the pions massive. To treat 
the theory in strong coupling we discretrize it by putting it on a spatial lattice with 
a lattice spacing of a = 1 fermi (/m) . R&WhH studied the physically relevant case of 
strong coupling: A = 20, f = 87.AMeV, and h = (119MeK)^. With these values the 
pion decay constant = 92.5MeV, the mass of the pion is = 135MeV, while 
that of the sigma is = GOOMeV. 

(The value of m„ = QOOMeV is traditional is the two flavor model. For the 
realistic case of three flavors in the next subsection, we shall instead identify the 
sigma field with a heavier particle, near IGeV in mass. For now I just note that a 
heavier sigma just makes the coupling even stronger than A = 20.) 

To understand the dynamics better, we also considered the model at extremely 
weak coupling, A = 10^^. In working at weak coupling we tried to keep at least some 
of the physics constant by holding the lattice spacing a and the pion decay constant 
fixed. This implies that the pion and sigma fields become light in weak coupling: 
for A = 10'^ and U = 92.5MeV, we obtain v = ST.iMeV, h = (2.03MeV^)^ 



m-TT = .3MeV and ma = l.SMeV. In other words, weak couphng means that the 
potential is fiat. 

We assume that the evolution of the system is "quenched". The idea is that 
the system cools so rapidly that the system finds itself in a state typical of high 
temperature, even though the evolution forward in time is by means of the equations 
of motion at zero temperature. We do not address the question of whether or not 
this is realistic for heavy ion collisions; it certainly is the extreme limit of a range of 
possibilities. We begin with a high temperature state which is chirally symmetric, 
($) = 0; to simulate the effect of fluctuations in the high temperature initial state, 
we distribute the fields as gaussian random variables with ($) = 0, ($^) = and 
($2^ = v'^/lfm^ following R&W.0 Pion domains can form in a quench because the 
chirally symmetric initial state, ($) = (0,0), is unstable against small fluctuations 
in the T = potential.il@ In essence, the system "rolls down" from the unstable 
local maximum of V^($) towards the nearly stable values $^ = -Lr^ ■ This process 
is known in condensed matter physics as spinoidal decomposition.El Long-lived DCC 
field configurations with tt 7^ can develop during the roll-down period. The field will 
eventually settle into stable oscillations about the unique /i 7^ vacuum, ($) = (/^r, 0). 
Oscillations continue until interactions eventually damp the motion through pion 
radiation. 

In the the Hartree approximation the equations of motion for the Fourier compo- 
nents of the pion field areilBS 

^vf, = (A(.^ - ($^)) - e) vf, . (2) 

Even in the Hartree approximation we see that field configurations with ($) = 0, 
(•l*^) < f ^ and momentum k < \fXv are unstable, and grow exponentially. Conversely, 
modes with larger momentum, k < \/~Xv are stable, and do not grow. Of all of the 
unstable modes, the constant mode, with = 0, grows the fastest; its growth has a 
natural time scale, which is 

Tsp = {X{v^ - {e))}-"^ ~ V2/m.. (3) 

In strong coupling, A = 20, this time is .'bjmjc. The exponential growth of the 
unstable modes continues until ($^) reaches w ^, when it begins to oscillate about the 
stable vacuum. Rajagopal and Wilczek found that the power oc vr'*^7r'*_^ in the low 
momentum pion modes indeed grows when the exact classical equations of motion 
are integrated for A = 20. Following Boyanovsky et aH, one can estimate how 
long it takes domains to grow; not suprisingly, this time scale is on the order of the 
spinoidal times, r^p. After this time, the fields sense that there is a stable minimum to 
the potential, and so from the nonlinearities in the potential, the growth of domains 
shuts off. Consequently, if the coupling were weak, so the sigma meson is light, then it 
takes long time to roll down to the bottom of the potential, and domains have plenty 
of time to grow. For the realistic case of strong coupling, however, the roUdown is 



very rapid, and so the domains are small, ~ Ifm. This is what one would expect from 
a system in which the Compton wavelength of the pion at rest is of order ~ Ifm. 

We confirmed this intuition through numerical simulations. An asymmetric lat- 
tice geometry of dimensions 10^ x 40 was chosen in order to study the issue of domain 
size in the one, long direction; we introduce the average of the pion field over the 
transverse dimensions, TiL{z,t) = X^x.y ^(^5 Z/5 ^5 ^)/(10)^- Initially, vr^ starts out com- 
pletely random. In weak coupling, by times of 100/m domains — regions in which 
the field is slowly varying about some nonzero value — are evident. In contrast, in 
strong coupling by times of 100/m the pion field is always quite random, oscillating 
with small amplitude about zero. We found that the pion correlation function is long 
ranged in weak coupling, but short ranged in strong coupling.l§n2l 

The most important quantity is the distribution in 7^3. We assume that we 
can extract this ratio from the classical pion fields in the most naive way possible, 
by computing 7^3 = ((vr3)^)/(((vr^)^) + ((vr^)^) + ((Tr^)^)). Histograms of 7^3 were 
obtained by evolving 200 independent configurations forward in time to t = 150 fm 
in weak coupling, and times of t = 30 fm in strong coupling. In weak coupling the 
distribution is far from binomial, and is approaching the DCC value of 1/(2^^7^3). 
In strong coupling, however, the distribution is clearly binomial, peaked about the 
expected value of 1/3. 

We have performed numerical simulations at other couplings to see where the 
crossover from a binomial to a DCC distribution occurs.^ For our lattice of 10^ x 40, 
it appears that a DCC distribution only arises at very weak coupling, A ~ 10~^, as 
expected from the domain size estimate of (2). 

In summary, in a realistic two flavor model we find that following a quench, DCC's 
are only produced in small, pion sized domains. There appears to be one way for large 
DCC's to be produced, and that is if there is a light particle about, whose Compton 
wavelength automatically provides a size for the DCC. But a quench alone does not 
seem to give us a large distance scale. 

3. 2. Three flavors 

In the next subsections I analyze the phase diagram of QCD^ This is similar to 
the analysis of sec. 2, except that I work not with 2 but the realistic case of 2 + 1 
flavors. The price I pay is that the contribution of the vector fields is dropped. Again, 
implicitly I speak of nonzero temperature; what happens at nonzero baryon density 
(if the temperature is small) may well be very different. 

From the previous discussion, remember that for three colors, the order of the 
phase transition appears to depend crucially on the values of the quark masses. This 
is definitely special to three colors. For example, consdier the limit of infinitely many 
colors, Nc 00. The basic physics can be understood simply by remembering that 
there are of order ~ gluons versus order ~ A^^^ quarks. Now assume that the 



deconfining transition is of first order aX Nc = oo; now while I certainly cannot prove 
it, the most natural possibility is that the first order deconfining transition alwaYS 
dominates the chiral transition, so that both take place at the same temperature.l2j 
This need not be the case, for it is logically possible for the chiral transition to occur 
at a higher temperature than the deconfining transition. (The opposite possibility, 
of a deconfining temperature which is higher than the chiral transition, probably 
contradicts general theorems on realizations of chiral symmetry, but as of yet no 
strict proof has been given.) Nevertheless, as the deconfining transition the free 
energy goes from being of order one, from hadrons, to of order N^, from gluons, and 
it is difficult to imagine that this in and of itself doesn't trigger the chiral transition 
at the same time. 

Thus we see that at least in this instance, three colors is not close infinity:0ii 
for three colors and "2 + 1" fiavors, while the transition appears to be of first order 
for the deconfining transition without dynamical quarks, Ad = oo and for the chiral 
limit, = 0, these lines do not meet — there is a gap, with QCD somewhere in 
between. 

Lines of first order transitions typically end in critical points, so it is natural to 
ask about the two critical points. Consider first working down from infinite quark 
mass. As Ai decreases from = cx), the line of deconfining first order transitions 
ends in a deconfining critical point. By an analysis similar to that in sec. (3.2) below, 
one can show that at the deconfining critical point, correlation functions between 
Polyakov lines are infinite ranged, and that it lies in the universality class of the Ising 
model, or a Z2 spin system, in three dimensions. 

The other limit is to go up from zero quark mass. For three fiavors the chiral 
phase transition is of first order at m = 0, so as m increases, the line of first order 
transitions can end in a chiral critical point. In the following we show that for 2 + 1 
fiavors, at the chiral critical point the only massless field is the sigma meson. Thus 
again the chiral critical point lies in the universality class of the Ising model, or a 
Z{2) spin system. 

We stress, however, very different fields becomes massless at each critical point; 
the two critical points are not naturally related to each other. At least for three colors 
this seems to be the there is a large gap in between the two critical points, in 

which there is no true phase transition, only a smooth crossover. Nevertheless, our 
results indicate that even in the region of smooth crossover — which includes QCD 
— that very interesting and unexpected physics can be going on. 

Because we are interested in making contact with QCD, we begin in the next 
subsection by fitting the scalar and pseudoscalar mass spectrum in QCD to the results 
found in a linear sigma model. This is an old story^ and relatively straightforward. 
The classification of the chiral critical point then follows directly. When we then 
try to connect the two, however, we find a surprise. Numerical simulations of lattice 
gauge theoryc3 indicate that as a function of the mass parameter A4, QCD is not far 
from the chiral critical point. We show that if QCD is in fact near the chiral critical 



point, then that, and the nature of the fit to the zero temperature spectrum, makes 
very interesting predictions about the behavior of QCD at nonzero temperature. 

3.2.1. The QCD spectrum at zero temperature 

In this subsection we use a hnear sigma model to model the full theory with three 
quark flavors.Silil0'il We ignore the contribution of vector mesons, which will be 
incorporated at_a later date.ta This type of analysis was first carried out by Chan 
and Haymaker;E3 re£ent analyses, with similar results, are given by Meyer-Ortmanns, 
Pirner, and Patkos,E3 and by Metzger, Meyer-Ortmanns, and Pirner.Ea For three quark 
flavors the sigma field $ is a complex valued, three by three matrix, proportional to 
the quark fields as $ ~ '(lieft bright- In flavor space it can be decomposed as 

8 

= E K + ^^«)^"- (4) 

a=0 

(ti . s are the generators of the S't/(3) algebra in the fundamental representation, while 
to is proportional to the unit matrix. Normalizing the generators as tr{tatb) = 6°'''/2, 
to = l/\/6.) The fields cTq are components of a scalar {J^ = 0+) nonet, those of TCa 
a pseudoscalar (J^ = 0~) nonet. The latter are familiar: 7ri_2,3 are the three pions, 
denoted as vr without subscript; the 7^4,5,6,7 are the four kaons, the i^'s, while vrg and 
TTo mix to form the mass eigenstates of the rj and t]' mesons. Following our notations 
for two flavors, we define the components of the scalar nonet analogously: we refer to 
c"i,2,3 as the cr,r's, to as the ai^-'s, while ag and ctq mix to form the and cr^/. 

What is the correct identification of the a meson for three flavors will be one of the 
principal questions in the following. 

This multiplicity of eighteen fields is to be contrasted with the usual two flavor 
sigma model considered in the previous subsection, which had only three vr's and one 
a meson. The increase in the number of fields is because we are allowing for ?7a(l) 
rotations; to describe two flavors in this way requires eight fields: the three vr's, a,,, 
the three cTtt's, and the rj. For rtf flavors, then, in general 2n^j fields are required to 
describe the lowest lying scalar and pseudoscalar multiplets. 

The effective lagrangian for the chiral $ field is taken to be 

C = tr - tr (h{^ + <l>^)) + /i^ tr (<l>"f$) (5) 

-\fQc(det{^) + rfet($^)) + {gi - (tr^^ + ?> g2 tr (^^^ 

This is very much like the lagrangian for two flavors in (1), except that now new 
couplings are possible. This lagrangian involves six parameters: two for the back- 
ground field H (one for the up equals the down quark mass, one for the strange quark 
mass), the mass parameter /i^, an "instanton" coupling constant c, and two quartic 
couplings, gi and g2. At large $ the potential is bounded from below if the quartic 



couplings satisfy gi > and gi + 2g2 > 0. (The bounds on the quartic couphngs can 
be understood by taking two extreme hmits for $: the first is $ proportional to the 
unit matrix, the second is where $ only has one diagonal element which is nonzero.) 

To understand the symmetries of this lagrangian, consider the theory becomes 
less symetric as various couplings increase from zero. When H = c = g2 = 0, the 
theory only involves the quadratic invariant tr($"l'$) = J2a=o (^a + ^a)> and so has 
the very large symmetry group of 0(18). When g2 7^ 0, the symmetry reduces to 
SU{3) X SU{3) X f/(l), with the U{1) that for axial fermion number. The effects 
of instantons, or more generally topological fluctuations and the Adler-Bell-Jackiw 
anomaly, generate a nonzero value for the instanton coupling constant c, and reduces 
the symmetry to SU{3) x SU{3). Lastly, spontaneous symmetry breaking, where 
($) ~ 1, occurs if the mass squared, /x^ is negative; actually, due to the determinental 
term, spontaneous symmetry breaking can even occur for some values of positive fi"^. 
Since ($) ~ 1, the symmetry is then reduced to the usual SU{3) symmetry of isospin 
plus strangness. 

For the background field H we take H = Hq tQ — ^/2hs tg- Working out the algebra, 
ho and are proportional to the current quark masses as 

Mup = Mdown ~ ho - hs , Mstrange ~ /^O + 2 /ig • (6) 

With the octet magnetic field hg ^ 0, then, only remaining symmetry is merely the 
usual SU{2) symmetry of isospin. Effects from isospin breaking, as from Aiup 7^ 
■Mdown, are negligible for our purposes, since every other mass scale is so much larger. 
(Not to mention the inclusion of electromagnetic effects.) Note that H has dimensions 
of mass cubed; hence two powers of some fixed QCD scale enter to make up the mass 
dimensions in (|^). I assume that these two powers of a QCD scale are truly fixed, and 
do not vary with, say, temperature. This appears to give a behavior with nonzero 
temperature which is consistent with other analyses, but I must confess that I do 
not not know how to derive the value of the background field from the fundamental 
theory in a way in which these two powers of the QCD scale would be manifest. n 

The analysis of the effective lagrangian is uniformly at the simplest level of mean 
field theory. I assume that in the true vacuum there are nonzero expectation values 
for CTg and o"o, 

c^o So + ctq , 0-8 + OS , (7) 

and then expand in powers of cXa and it a- To linear order there are two equations of 
motions, which fixes the values the two vacuum expectation values, Sq and Eg, as 

ho = Eo - c (e^ - S^) + 4g2T.l(j:o + l^8) , 

hs = Eg + cSg (So - Sg) + 2g2 Sg (Sq + Sg) (^Sq + ^ Sg) . (8) 
A second mass parameter, M^, has been introduced; it is related to the /x^ of (|^) as 

= + + 2S2) . (9) 



Expansion of the effective lagrangian to quadratic order gives the masses of the pseu- 
doscalar octet, 

ml = M2 - c(So + 2S8) - 2(72^8 f So + ^ S ^ ^" " ^ 



^8 



mi, = - c (So - Sg) + Sg (So + 5S 



2 / So — Sg 



So + Sg/2 
2 



m^^ = - c (So - 2Sg) + 2(72 Sg (^So + ^ Sg 

"^^sTTo = Sg (^c - 2^2 + ^ Sg 

m^^ = + 2cSo . (10) 

Notice that m^^, '^T'^gTro' ^^"^ "^^o ^^fer to the octet and singlet parts of the pseudoscalar 
nonet, which mix to give the masses of the rj and m,, and m^/. In the chiral hmit, 
ho = hs = 0, the chiral condensate is SU{3) symmetric, Sq 7^ and Sg = 0. In this 
limit it is easy to see that the tt, K, and 77 mesons are all massless; the r/' meson 
remains massive because of the instanton coupling c. 

For the scalar nonet, 

ml^ = + c (So + 2Sg) + 2(72 f S^ - SSoSg + i Sg 



2 

mi^. = M' + c(So-S8) + 2(72 (So + Sg) ( So + ^ S, 
ml^ = /x^ _ c (So - 2S8) + (71 (S^ + eS^) + 2(72 ( S^ + SSoSg + ^ S^ 



2 - -v^Sg (c +2(7iSo + (72(4So + 3Sg)) 



= _ 2cSo + (71 (3S2 + 2S2) + Ag^T^l . (11) 



For the scalars, m^^, "^ctso-o' ^'^'-^ "^^o "^^^ Si^e mg-^ and m^^,. All particles in the 
scalar nonet are expected to remain massive even in the chiral limit. 

We also need the pion and kaon decay constants, 

U = \ l (So-Sg) , fK = \ l fSo + ^Sg) . (12) 



3 ' , ... V 3 V 2 

The expressions in (p!0|)-([T^) agree with previous results, up to differences 

in normalization of the coupling constants and vacuum expectation values in (^ and 
(1^). There is one unexpected detail: naively one would expect that the mass and the 
quartic coupling gi would enter independently. In fact, for the equations of motion 
in (|^), the masses of the entire pseudoscalar nonet in (|10|), and the masses of half the 




scalar nonet — for the and the ax in (|TT]) — these two parameters only enter in 
a certain combination, through the single parameter of (^. This means that we 
can fit to the pseudoscalar spectrum, and so fix M^, and yet still be free to vary /i^ 
(or conversely gi): the only change will be to alter the masses of the an and the a^i. 
This technical detail will play an important role in what follows; although there must 
be some simple group theoretic reason for it, I don't know what it is. 

There is some freedom in deciding how to fit the parameters of the linear sigma 
model. Various kinds of fits are given by Mever-Ortmanns, Pirner, and Patkos,E2l and 
by Metzger. Meyer-Ortmanns, and Pirner.S Following the experience of Chan and 
Haymakerc2l we do not fit to the entire pseudoscalar mass spectrum for the vr, t], 
and 7]' mesons, since it turns out that the kaon mass is fairly insensitive to the ratio 
of vacuum expectation values, Eg/So- 

Because of this, we leave the ratio Sg/So as a free parameter, and fit just to the 
pion decay constant and to to the masses for the tt, t], and rj' mesons. I assume 
the values: 



= 93MeV , = 137 MeV , = 547 Mel/ , m^, = 95SMeV. (13) 
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Figure 7: Coupling for the anomaly, c, 



vacuum expectation values, s 
to s. 



in MeV, versus the ratio of octet to singlet 
Sg/So- This coupling constant is rather insensitive 



Then one can examine the sensitivity of various quantities to s = Sg/So- This 
is easiest of understand graphically: in fig. (6) I show how the kaon mass changes 
with s: for small s it is a little high, but not terribly so. In fig. (7) I show how the 
coupling for the anomaly changes with s; again, there is some variation but nothing 
dramatic. 

In fig. (8) I show how the coupling constant g2 of the linear sigma model changes 
with s; we see that for large s ~ .3 the coupling g2 becomes rather small. Analysis 
shows that while g2 becomes negative for s > .3, its value is very small, < 1., so that 
for any reasonable values of gi, the potential remains bounded from below. Thus 
values of s < 1. cannot be excluded on such grounds. 

There are two quantities, however, which are very sensitive to the value of the 
octet to singlet vacuum expectation values, s. The first is the kaon decay constant, 
fx- This is illustrated in fig. (9), where we see that small values of s are preferred: 
the experimental value is fx = 113 MeV".E2l Since fx isn't too far from given the 
definition of fx in (IHD, it is obvious why a small ratio of Eg/Eo is preferred. The 
second quantity is the mixing angle between the rj and the rj', Or^r]'- Experimentally 
the situation is a bit unclear: radiative decays favor a value of 6*^^/ = — 20°, while fits 
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Figure 8: The coupling g2 of the hnear sigma model versus the ratio of octet to singlet 
vacuum expectation values, s = Eg/So. While g2 becomes much smaller for large s, 
this is rather innocuous. 



to the mass spectrum favor ^^^/ = — 10". From fig. (10), however, we see that this 
uncertainty doesn't really affect our fit, for either value favors low values of s. From 
fig. (9) and (10), we find that a good fit is obtained for s = .1. The parameters of 
this solution are 



+ {U2MeVf , c = 1920 Me\/ , §2 = 30. 



127 MeV , Sg = 13MeV , ho = {290 MeV)^ , hs = {281 MeVY 



For this fit with s = Sg/Sg ~ .1, while 

rriK = 516 MeV , fx = 109 MeV , 9r,n' 



10.4° . 



(14) 



(15) 



The kaon mass is a bit high, rriK = 516 MeV instead of the (average) experimental 
value of 497 MeV". 

For the values of (0), from (|]) the ratio of the strange to up (= down) quark 
masses is 

■Mstrange _ (^0 + 2/lg) 



M 



up 



{ho - hs) 



32 . 



(16) 



. 05 



. 15 



0.2 



0.25 



Figure 9: Kaon structure constant, fx, in MeV, versus the ratio of octet to singlet 
vacuum expectation values, s = Sg/So- Small values of s are favored, s ~ .1. 

This is larger than the often quoted value of ~ 20 because we have explicitly allowed 
for a vacuum expectation value that is not SU (3) symmetric. Eg ^ 0. 

If one fits to the kaon mass, then Sg/So is constrained to be = .19, with c = 
1790 MeK and g2 = 4.98; the results for fx and 9r,n' are worse, fx = 127 MeV, 
while Oyri' = + 1-°- The freedom to play with the ratio Sg/So and its effects will be 
discussed elsewhere. 

The values in (0) give unique predictions for two masses in the scalar nonet: 

m^^ = 1177 MeV" , = 1322MeV. (17) 

There are observed stateJl^ with these quantum numbers: the ao(980) and the i^'g (1430), 
respectively. 

For the purposes of illustration, in fig. (11) I illustrate how the masses of the cTtt 
and the ax change with s. For small values of s, the ax is heavier than the cTtt, as 
naive intuition would suggest. What is amusing is that for s > .2, the ax, which 
after all is strange, becomes lighter than the 0",^. For other reasons, I do not favor 
such large values of s, but it is interesting that such a flip of mass values can occur. 

The identification of the with the cr,r is not obvious (see pg. VII. 21 of ref. |[40|| ). 
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Figure 10: The mixing angle between the rj and the r)', O^jn', versus the ratio of octet to 
singlet vacuum expectation values, s = Eg/So- Small values of s are favored, s < .1. 



It is reasonable to expect that the width of the o"^ should be large, while the Oq is 
observed to have a small width ~ 57MeV. (The Kq is a broad resonance, with a 
width ~ 287 Me V".) While this may result from a large mixing with other channels, 
such as alternately it is possible is that the ao(980) is not what we term the 

a-jr, but a kind of KK mplecule.E^ If true, presumably the state which we term the 
cTtt is heavier than 1 Gd^-cUcS 

There is no unique prediction for the masses of the remaining members of the 
scalar nonet, the o"^ and the ar^'. As remarked, the masses of all other fields only 
depends upon the parameter in (j^). Thus at constant we can vary gi and 
change only rricr^ and rricr^,. This is illustrated in fig. (12). 

Like the rj and the t]' mesons, the and the a^/ mix; the mixing angles are 
usually of order 30''. This is close to ideal mixing, where the a^/ is primarily uu + dd, 
while the cx^ is primarily ss. This mixing helps to understand the differences in the 
mass spectrum between the scalars and pseudoscalars. For the pseudoscalars the r] 
is primarily an SU{3) octet, and the r]' primarily an SU{3) singlet; the singlet r]' is 
then much heavier than the octet f] because of the anomaly. Indeed, the r] and r]' are 
primarily SU (3) because of the anomaly. For the scalar particles, instead of being 
SU (3) eigenstates, the cr^ and the cr^/ are much better described as flavor eigenstates. 



Figure 11: The masses of the scalar pion, a^^, and the scalar kaon, ax, versus the 
ratio of octet to singlet vacuum expectation values, s = Sg/So. 



Since the a,,' is (largely) made up of up and down quarks, then, and the cr^ is (largely) 
made up of strange quarks, it is natural that the for the scalars, the is heavier 
than the cx^'. 

(I note that in strong coupling, gi oo, m„^, approaches a constant value, while 
"^o-^ ~ y/9i ~^ C)0. This is clearly a special property of taking the up and down quark 
masses different from the strange quark mass.) 

I identify the cr^' and the cr^ with the observed state£^ with the same quantum 



numbers: the /o(975) and the /o(1476), respectively. With the parameters of (|T^ 
from fig. (12) requiring mg-^, = 975MeV fixes 

m^^ = U76MeV , gi = 40. , /i^ = -{A92MeVf. (18) 

For these values, the mixing angle between the o"^/ and the o"^ is +28". The mass for 
the (T^ is a bit high, U76MeV instead of UOOMeV for the /q. 

Again, the identification of the cr^/ with the /o(975) is not evident (see pg. VII. 192 
of ref. 1^): the /o(975) has a width of only ~ 47Me\^, while the a^/ should decay 
into two pions with hearty abandon. (The width of the /o(1400) is large, 150 to 400 
MeV.) Either the narrow width is the result of a complicated multi-pole structure in 
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Figure 12: Masses of the cr^ and the cr^/, in MeV, versus the couphng gi of the hnear 
sigma model for three flavors. Notice that mo- , = for (71 ~ 4.; this will be used in 
sec. (3.2.3). 

the KK channel,^ or the /o(975) is not the a,,/, but a KK molecule.E Indeed, since 
in this model the a^/ is primarily composed of up and down quarks, and very little 
from strange quarks, it is impossible to understand why it should decay into KK in 
the first place. As shall become clear from our later arguments, for our purposes all 
we really need to know that the a,,/ is not light; i.e., ttLct^^, > 1 GeV. 

So where does the a meson at GOOMeV, as we assumed in sec. II, come from? 
This is derived from models of nucleon-nucleon forces, assuming that the force is 
provided by single particle exchange. An isosinglet 0"*" particle, the a, is needed at 
GOOMeV to provide sufficient attraction. It is reasonable to view this a simply as 
an approximation to two pion exchange. Notably, there is no such resonance seen in 
the phase shifts of vr — tt scattering below IGeV, pg. VII. 37 of ref. |^0|. For the 
sake of argument, if we do assume that rrif^ , = QOOMeV, we find that the coupling 



gi decreases, to gi = 16.2, whilst m^^ = U12MeV and /x^ = +{384:MeVf. 

Other fits to the a^/ have also been given. Bertsch et al0 find that data on kaon 
decays and ttN ttttN scattering can be well fit with m^^, ~ SOOMeV. 

1 now show that understanding where the a^' lies is not merely to interest to 
afficianadoes of the linear sigma model: it is crucial to understanding the nature of 
the phase diagram, as a function of the current quark mass Ad. 



3.2.2. The chiral critical point 



I now turn to the nature of the chiral critical point; the above classification of 
the zero temperature spectrum is not required in this subsection, but will be used 
in the next. I deal exclusively with mean field theory, where the effects of nonzero 
temperature are incorporated simply by varying the mass parameter /x^. This is 
correct in the limit of very high temperature, but should be qualitatively correct at 
all temperatures. 

I begin with the SU{3) symmetric case when = 0. For a constant field Eq, the 
la grangian reduces to the potential for Sg, 

C = -hoJ^o + ^f^'^l - + (19) 

This model has precisely the same phase diagram as that for the transition between 
a liquid and a gas. For zero background field, = 0, the instanton coupling "c" is a 
cubic invariant, and so drives the transition first order. As ho increases the transition 
becomes more weakly first order, until at h^ = /ig^** the line of first order transitions 
ends in a critical point. For ho > h^"^ there is no true phase transition, just a smooth 
crossover. 

There is a critical point where the potential in Sq — Sq^** is purely quartic, 

i-icrit — gi I ^ I • l^uj 

This is determined algebraically, by requiring that about the critical point, in Sq — 
Eq''**, the first, second, and third derivatives of the potential vanish. This is three 
conditions, which fixes the values of the magnetic field, the vacuum expectation value 
of the So field, and the mass parameter, 

3 2 

h'^rit _ _^ ycrit _ 2 _ (r),^ 

The mass spectrum at the critical point is easy to work out by going back to the 
original action and recomputing: 

2 2 '- 2 1 n 2 

= niK = = y , m^, = lOm^ , 



ml^ = ml^ = ml^ = (7+18^) ml , = Q. (22) 



9i, 

Of course the ar,' is massless because we have dehberately tuned ourselves to sit a 
the chiral critical point. As we work for now in the limit of SU{?)) symmetry, there 
are two octets of pseudoscalars and scalar fields. It certainly is extraordinary to find 
even an isolated point where the a^/ is lighter than any other field, even the pions! 

Because only one field, the o"^/, is massless at the chiral critical point, the similarity 
to the liquid gas phase transition extends to the universality class: it is that of a Z{2) 
spin system, or the Ising model in three space dimensions. Thus the chiral critical 
point is precisely analogous to critical opalescence in a liquid gas system. 

The classification of the universality class still applies when the background mag- 
netic field is not SU(?>) symmetric, ^ 0. In this case I couldn't solve the equations 
analytically, and had to resort to numerical methods (which are, however, elemen- 
tary). The difficulty is simply that three equations must be set to zero to find the 
chiral critical point, for the first, second, and third derivatives of the potential, in the 
direction of the massless mode. (And of course for the first derivative of the potential, 
in the other direction.) By solving ( p!oD and (|l3); I find that at least for parameters 
about those of the previous section, there remains a single, massless field at the chiral 
critical point, the cr^/, with the universality class that of the Ising model. Of course 
for /ig 7^ the o"^/ field does not remain a pure S'?7(3) singlet, but mixes to become 
part octet. 

I expect this conclusion to hold for arbitrary values of the quark masses (back- 
ground field). Indeed, I conjecture that in general, for any critical end point — be 
it chiral, deconfining, or whatever — that the universality class is always that of 
the Ising model, which is 2'(2). I have no profound reason for this, but simply on 
algebraic grounds, at best one can get the appropriate derivatives of the potential to 
vanish along one direction, but not in more than one direction; if one tried in more 
than one direction, there would be more conditions on the potential than variables 
to satisfy them. 

Indeed, although I did not comment on it at the time, the identification of the 
Gr^i as the massless mode of the chiral critical point could have been anticipated from 
the results of the previous section. Looking at fig. (12), we see that for very small 
values of (^i, (?! ~ 4., that the mass of the cr^/ goes to zero. At the time, it seemed 
merely a peculiarity of the fit, of no greater significance. And yet it shows that even 
at zero temperature, it is possible to obtain a massless a^'. 

Assuming, then, that only the o"^' is massless at the chiral critical point, I can 
flesh out the entire phase diagram as a function of M.up = M-dmun versus M-strange-, as 
proposed in fig. (1) of ref. 0. The basic idea is that along the entire line of chiral 
critical points, only the cr^/ is massless, but that the singlet-octet ratio in the cr^/ 
changes. I just showed that at the ^[/(S) symmetric point, M.up = -Mdown = -Mstrange, 
the ar,' is pure singlet. Decreasing strange to the critical point where strange = 
and A4up = M.down 7^ 0, rather obviously the a,,/ must become entirely strange. 



a„i ~ ss. 



The opposite limit, Aiup = ■M.^own = 0, is more familiar. Assume that the chiral 
phase transition with two massless flavors is of second order?^ the universality class 
is that of an 0(4) critical point. The chiral phase transition will remain of second 
order as Mstrange decreases from Aistrange = oo. This line of second order transitions 
must end in a special value of Mstrange = A^stmnge! fo^^ Mstrange < MZl^g^, the chiral 
transition is like that of three light flavors, and so first order. Wilczeke^ observed that 
at Aistrange = ■At'^^^^g^, the chiral transition is in the universality class of an 0(4) 
ih-critical point. I would say that when A4up = -Mdown = and Aistrange = -^strange 
the (j^/ is a pure SU (2) state, with ar,' ~ uu + dd; further, because Aiup = Aidown = 0, 
the pions are also massless, so the universality class of this point is not that of a Z[2) 
critical point, but an 0(4) tri-critical point. This contradicts my previous argument 
(that only the cr^' is massless at a chiral critical point) because when Aiup = Aidown = 
0, that is itself the endpoint of a line of chiral critical points, and thereby special. 

The analysis can be further extended to the case where there is a gap between the 
first order transitions for the deconfining and chiral transitions, but where the chiral 



transition is of first order for two, massless fiavors. Then fig. (1) of ref. would 
be modified, with a band of first order transitions about the axis Aiup = Aidown = 0. 
These first order transitions would end in chiral critical points, in the universality 
class of Z{2) from the presence of massless a^/ fields. However, the latest lattice data 
appears to indicate that this is not the case: that for two fiavors, the chiral transition 
is of second order. 



3.2.3. The chiral critical point and QCD 

To model the effects of nonzero temperature, in mean field theory it is only 
necessary to increase the mass parameter /x^ of (|^). Doing so, the resulting spectrum 
— computed from (^-([11]) — is rather unremarkable, and look very similar to the 
graph in fig. (1), except that now there are many more fields about; this is illustrated 
in fig. (5). Remember fig. (1) refers to two fiavors, including vector mesons, while 
fig. (5) is for three fiavors, without vector mesons; also, in fig. (1) I took mo-(O) = 
rricr^ = 600 MeV, while in fig. (5) mo-^, = 1 GeV. Yet the qualitative behavior of the 
states is the same: the only two states which change significantly with increasing /i^ 
(which corresponds to increasing temperature T) are the pions and the cx^/. For three 
fiavors, in fig. (5), the cx^/ comes down in mass from ~ 1 GeV, while the pion moves 
up in mass, until they meet at about ~ 400MeV^; then both increase together. Since 
the cr,,' is coming down in mass from 1 GeV, it is natural to find that it meets up 
with the TT at a much higher mass, ~ 400 Mel^, than before; if the a starts out at 
600 MeV, then it matches up with the vr at ~ 200 MeV. 

The limitations of mean field theory are clear. Including the effects of nonzero 
temperature by only changing /x^ corresponds to the addition of an SU{3) symmetric 
mass term ~ T^tr($^$). While this type of term is dominant at high temperature. 



at low temperature things are surely much more complicated. Calculating just the 
thermal effects to one loop order for the model of (^) would be an arduous task. 
Nevertheless, we observe that mean field theory is in qualitative agreement with 
more sophisticated analyses in effective sigma models, instanton models, Nambu Jona- 
Lasino models, etc.. With one exception, these models a// come to the same conclusion 
as mean field theory: in the chiial limit, the chiral transition is of second order for two 
flavors,cl'E3 first order for three,EZl while for the nonzero quark masses of QCD there is 
only a smooth crossover .l3caE3 I stress that this agreement is rather remarkable, since 
very different methods of calculation are being employed. 



The one exception is the work of ref. who find a second order chiral phase 



transition for three massless fiavors. These works, however, use an approximate 
solution to the Schwinger-Dyson equations. This approach does not include the axial 
anomaly, and so in this approximation should give a second order phase transition. 
Thus the one discrepancy is understandable on general grounds. 

We can now use the fit to the zero temperature spectrum to ask a more detailed 
question: how far is QCD from the chiral critical point? The data of ||3^ indicates 



that as a function oi Ai = M.up = -M-down = M. strange-, QCD is only about a factor 
of two from the chiral critical point. 

Let us assume that the chiral phase transition is of first order for three, massless 
flavors because of the presence of the instanton coupling ~ det{^). Then we assert 
that it is very difficult to understand why QCD is (relatively) so close to the chiral 
critical point. Varying the current quark masses is equivalent to varying the back- 
ground fields /iq aiid /ig. We then take (P)-(PT^, and search numerically for the values 
of ho = /iq'"** and hg = hg"^ which give a chiral critical point. In doing so we require 
that the ratio of up to strange quark masses is fixed; from (^, this implies that the 
following ratio of ho and hs is fixed to have the same value as at zero temperature: 

(ii±|i). 32. (23) 

{ho - hs) 

By varying ho, hg, /i^, Sq, and Eg, and otherwise using the values in (|14D and (pTsD , 



I determine where the critical point is. There are four equations which must vanish: 
two equations of motion, one equation to ensure the a^/ is massless, and lastly, one 
equation so that the cubic variation of the potential, in the direction of the Ur,', 
vanishes. It is this last equation that makes the task of determining the location of 
the critical point so difficult analytically. Numerically things are easy, using these 
four conditions and p3|) to determine 

h^"' = {62MeVf , h^"' = {60.4MeVf , fil^ = {183 MeVf , 

So = 14.5MeF , Sg = 2.7MeF. (24) 

I then compute the ratio between the current quark masses at the chiral critical point 
to those in QCD: 

= —u ^ = • 25 

Mup ho - hs 



Now we easily confess that this ratio, computed in mean field theory, is crude at best. 
Even so, numerical simulations^^ find that the ratio in ( p5D is not 1%, but 50% — 
mean field theory is off by almost two orders of magnitude! 

For the purposes of discussion, the ratio in (^) is not significantly larger for a 
lighter o meson. If we take m„ , = GOOMeV, one finds that A^^pY-^up = -06 instead 
of .01. This is still a long way from the value of ~ .5 found by the Columbia group.il 

Also, I note that Meyer-Ortmanns, Pirner, and Patkos0 had previously argued 
that (in our language) QCD is far from the chiral critical point. They did not make 
a precise estimate of this distance, though, as in (|25|) . 

This leads us to discard our initial assumption. In the chiral limit of three, 
massless flavors, there are two mechanisms for generating a first order transition. The 
first is the presence of the (cubic) instanton coupling, ~ (iet($).0ii From (^), this 
does not give a phase diagram that is even qualitatively correct. Therefore, perhaps 
the second mechanism is at work:§ in the chiral limit, the chiral jihase transition is of 
first order because it is a type of Coleman- Weinberg transition.caH 

In mean field theory quartic couplings are fixed and do not change with temper- 
ature. I use the term Coleman- Weinberg to describe a theory in which the quartic 
couplings change significantly with temperature (in condensed matter physics this is 
known as "fluctuation induced"). This phenomenon can be demonstrated rigorously 
in fourii'ill and 4 — eil dimensions; extrapolation to three dimensions, e = 1, may 
or may not be valid. Numerical simulations of SU{3) x ST/Xiilspin models in three 
dimensions suggest that they are Coleman- Weinberg like;E2E3E3 one then has to rec- 
oncile this with the expansion up from two dimensions, which predicts a second order 
phase transition.ii'ii 

(Incidentally, the claim of Chivukula, Golden, and Simmonsiil that the /3-functions 
in ref. (0) is incorrect is itself incorrect: the discrepancy is only apparent, from using 
a different normalization of the kinetic term. I thank Y. Shen for pointing this out.) 

We propose that the only way for the phase diagram presently observed in nu- 
merical simulations of 2 + 1 flavors to be true is if the quartic couplings of the effective 
linear sigma model run by large amounts. In 4 — e dimensions,E3o in the space of 
Qi and g2 the couplings tend to run most for large g2, from large to small values of 
gi for approximately constant g2. This type of flow is precisely what is necessary to 
obtain a massless cr^' at the chiral critical point. As I noted previously, this can be 
seen even from the parameters of the model at zero temperature, where a massless 
arji is obtained by letting (?! — > 4.. This hardly proves that such a speculative sce- 
nario is right; but at least couplings flow in the right direction, with features that are 
qualitatively correct. 

Of course, I confess that using primitive mean field theory is a gross caricature 
of reality. Nevertheless, the question of how far QCD is from the chiral critical point 
does not seem to have been asked previously. It is surprising that the question of 
exactly where and what the o"^' meson is has such a dramatic effect on the phase 



diagram of QCD. If nothing else, while are predictions are probably baseless, they 
are precisely testable in numerical simulations of lattice gauge theory. Hopefully we 
won't have to wait until the next millenium for simulations with dynamical fermions 
that are good enough to be believeable. 



3.2.4- Three flavors and DCC's 

Given the crudeness of present day numerical simulations for 2 + 1 flavors, I 
conclude with wild speculation: perhaps QCD is very close to the chiral critical 
point. A priori, there is no reason why QCD should be very close to the chiral 
critical point; there would have to be some magical conspiracy of mass scales. 

But if it is, then even if there is no true phase transition in QCD, there is a natural 
mechanism for the generation of a large distance scale at nonzero temperature, from 
a light arf'. As we alluded to in a previous work,Ea this large distance scale might 
be used to produce large domains of Disoriented Chiral Condenstates. This is not 
obvious, for even at the chiral critical point, only the arj' becomes massless, not the 
pions. The question is then: are large domains of DCC's generated when massive 
pions are coupled to a massless (or light) a^/? 

One way is by the type of isospin zero states constructed by Horn and Silveili 
and discussed recently by Kowalski and Taylor.ii These authors construct a state 
with total isospin zero for which distribution in isospin is broad, like that of a DCC. 
The question is then, is it reasonable for such states to be produced near the chiral 
critical point ?ii 
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